Міністерство освіти України
ІКТАМ
при Національному університеті "Львівська політехніка"
Кафедра : автоматики та телемеханіки
Курсова робота
з курсу
"Чисельні методи аналізу автоматичних систем"
на тему:
"Система статичного слідкування за швидкістю задаючого вала"
Тема 4, варіант 11
Львів 2003
Зміст
Завдання 3
Опис роботи схеми 4
Виведення системи диференціальних рівнянь 4
Метод Рунге-Кутта-Мерсона 5
Модифікований метод Ейлера 5
Реалізація автоматичного кроку 6
Блок-схеми алгоритмів. 7
Тексти програм 8
Результати 13
Графічні результати 5
Висновок 6
Література 6
Завдання
Двома заданими чисельними методами розрахувати з заданою точністю процес реакції заданої системи автоматичного керування (САК) на одиничний стрибок вхідного сигналу.
Розрахунок здійснити для свого варіанту САК, параметрів САК, параметрів перехідного процесу і параметрів процесу чисельного інтегрування.
Чисельні методи 1) з фіксованим кроком – Рунге-Кутта-Мерсона.
2) з автоматичним вибором кроку – Модифікований Ейлера.
Точність розрахунку не нижча 0.001, початкові умови – нульові.
1. Схема системи статичного слідкування за швидкістю задаючого вала.
2. Рівняння ланок
a) рівняння тахогенераторів
,
b) вимірювальна схема
c) електронний підсилювач
d) обмотка збудження ЕМП (електромашинного підсилювача)
e) короткозамкнута обмотка ЕМП
f) двигун
3. Завдання
Tm (сек)
0,006
КЕМП
20
T1 (сек)
0,2
КЕП
1.5
T2 (сек)
0,015
к1 (в∙с/рад)
1,01
С (рад/в∙с)
1
к2 (в∙с/рад)
1
Опис роботи схеми
На вхід електронного підсилювача подається різниця напруг U1 та U2 які є пропорційними величинам вхідного і вихідного сигналів відповідно. Сигнал, підсилений електронним, а потім електромашинним підсилювачами, подається на двигун, який обертає ротор тахогенератора ТГ2.
У рівноважному стані Wвх = Wвих, на вході електронного підсилювача нульова напруга. При появі , припиненні чи зміні Wвх на вході електронного підсилювача (ЕП) виникає напруга: Uвх=U1-U2 . Підсилена в Кu разів різниця напруг подається на обмотку збудження електромашинного підсилювача (ЕМП), який керує двигуном. Під дією цієї напруги в залежності від знаку зміни вхідного сигналу збільшується або зменшується швидкість обертання двигуна. Вал двигуна з'єднаний з тахогенератором, який генерує напругу, що подається на вимірювальну схему. Цей процес буде тривати доти , доки Uвх не стане рівне нулю , тобто Wвх = Wвих .
Виведення системи диференціальних рівнянь
З рівняння двигуна виражаємо Up :
і підставляємо його в рівняння короткозамкненої обмотки ЕМП, звідки знаходимо Uk:
.
Підставивши вираз для Uk, а також рівняння електронного підсилювача, вимірювальної схеми і тахогенераторів в рівняння обмотки збудження ЕМП, отримаємо:
Позначимо :
тоді .
Тепер, з отриманого лінійного диф. рівняння третього порядку можна скласти систему з трьох диф. рівнянь першого порядку:
Метод Рунге-Кутта-Мерсона
Метод Рунге-Кутта-Мерсона є модифікацією метода Рунге-Кутта
4-го порядку. Він реалізується за наступним алгоритмом:
За допомогою циклів з керуючою змінною j обчислюємо
Yj(n+1) = Yjn + ( K0j + 4K3j + K4j ) / 6,
де
K0j = h * Fj( Xn , Yjn )
K1j = h * Fj( Xn + h/3 , Yjn + K0j/3 )
K2j = h * Fj( Xn + h/3 , Yjn + K0j/6 + K1j/6 )
K3j = h * Fj( Xn + h/2 , Yjn + K0j/8 + 3K2j/8 )
K4j = h * Fj( Xn + h , Yjn + K0j/2 - 3K2j/2 + 2K3j )
Метод Рунге-Кутта-Мерсона також має похибку інтегрування h5 , але, на відміну від методу Руеге-Кутта 4-го порядку, для оцінки похикби відпадає необхідність робити подвійний перерахунок. Цей метод дає наближену оцінку похибки на кожному кроці інтегрування за формулою:
Rj(n+1) = ( -2K0j + 9K2j - 8K3j +K4j )/30
Модифікований метод Ейлера
Розв’язок в околі точки xm розкладається в ряд Тейлора
Якщо xm змінюється з кроком h, припустивши, що сітка x0,x1,...,xm є рівномірною, можна отримати наближене значення y(x) в точці xm+1...